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PREDICTION  OF  TIME-TO-GO  FOR  A  HOMING  MISSILE 
USING  BANG-BANG  CONTROL' 
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Department  of  Aerospace  Engineering  and  Engineering  Mechanics 
The  University  of  Texas  at  Austin 
Austin,  Texas 


Abstract 

The  flight  time  required  for  a  variable-speed  homing  mis¬ 
sile  to  intercept  a  zefo^ccelefatlon  target  in  two-dimensions 
is  determined  by  assuming  that  the  iiussile  normal  acceler¬ 
ation  is  bang-bang,  that  is,  maxintum  norih^_a^leration 
followed  by  zero  normal  acceleration  .where  the  swit^.Uiqe 
which  gives  intercept  is  to  be  determined.  For  those  cases 
where  intercept  does  not  occur  in  a  reasonable  time,  the 
flight  time  which  minimizes  the  miss  distance  is  used.  A 
tangential  acceleration  profile  is  assumed  for  the  missile, 
that  IS.  constant  positive  acceleration  when  thrusting,  con¬ 
stant  negative  acceleration  when  coasting,  and  a  given  en¬ 
gine  burn-out  time.  In  this  way,  the  velocity  of  the  missile 
becomes  a  known  function  of  time,  and  the  equations  of  mo¬ 
tion  can  be  solved  analytically.  Then,  an  algebraic  equation 
for  the  switch  time  for  intercept  or  the  final  time  for  closest 
approach  can  be  derived,  but  it  must  be  solved  numerically. 

The  time-to-go  algorithm  is  tested  in  a  six-d^ree-of- 
freedom  simulation  of  a  homing  missile  with  a  linear- 
quadratic  guidance  law  where  the  target  performs  two  maxi¬ 
mum  normal  acceleration  maneuvers.  At  each  sample  point, 
the  missile  velocity  vector  is  projected  onto  the  plane  of  the 
line-of-sight  vector  and  the  target  velocity  vector,  and  the 
time-to-go  is  calculated  for  a  planar  intercept.  This  time- 
to-go  is  used  to  calculate  the  gains  for  the  guidance  law. 
Results  show  that  this  time-to-go  algorithm  improves  the 
performance  (miss  distance)  of  the  missile  for  several  sce¬ 
narios  relative  to  the  range-over-closing-speed  sJgorithm. 

List  of  Symbols 

A  F\inction  defined  in  Eq.  (14) 

On  Missile  normal  acceleration  (ft/sec^) 

oi  Missile  tangential  acceleration  (ft/sec^) 

ai ,  Constants  defining  VX/  for  thrusting  flight 

02,  bj  Constants  defining  Vm  for  coasting  flight 

B  Function  defined  in  Eq.  (16) 

d  Miss  distance  (ft) 

i  Thrusting,  1;  =  1;  coasting,  k  =  2 

f  Time  (sec) 

t„  Initial  time  (sec) 

t;  Engine  cutoff  time  (sec) 

ty  Final  time  (sec) 
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t.  Switch  time  between  a„  =  Ori,u,  and  On  =  0  (sec) 

Missile  velocity  (ft/sec) 

Vr  Target  velocity  (ft/sec) 

X  Inertial  X  coordinate  (R) 

z  Relative  X  coordinate  (ft) 

Y  Inertial  Y  coordinate  (ft) 

y  Relative  Y  coordinate  (ft) 

ff  Missile  velocity  orientation  (deg) 

<t>  Target  velocity  orientation  (deg) 

Subscripts 

c  Engine  cutoff 

/  Final 

o  Initial 

p  Point 

s  Switch 

Introduction 

A  guidance  law  of  current  interest  for  bank-to-turn  hom¬ 
ing  missiles  is  the  linear-quadratic  guidance  law  which  con¬ 
tains  proportional  navigation  as  a  particultir  case  (see,  for 
example.  Ref.  1).  In  order  to  implement  this  guidance  law, 
an  algorithm  for  predicting  time-to-to  is  needed.  The  sim¬ 
plest  time-to-go  formula  is  range  divided  by  closing  speed 
and  is  valid  for  a  constant-velocity  missile  and  target  on  a 
collision  course.  This  formula  has  been  improved  in  Ref.  2 
by  accounting  for  the  missile  longitudinal  acceleration. 

Unfortunately,  the  linear-quadratic  guidance  law  tends 
to  drive  the  missile  and  the  target  into  a  homing  triangle 
in  which  range  and  closing  speed  become  unobservable,  In 
Ref.  1,  a  linear-quadratic  guidance  law  for  dual  control 
(intercept  and  estimation  enhancement)  has  been  proposed. 
This  guidance  rule  moves  the  missile  away  from  the  homing 
triangle  improving  estimation  but  making  the  time-to-go 
algorithm  invalid. 

The  purpose  of  this  study  is  to  develop  a  time-to-go  al¬ 
gorithm  which  is  valid  for  intercept  geometries  which  differ 
greatly  from  the  intercept  triangle,  such  as  encountered  with 
dual  control.  This  is  accomplished  by  assuming  tiiat  the 
missile  normal  acceleration  history  is  maximum  normal  ac¬ 
celeration  followed  by  zero  normal  acceleration  (bang-bang) 
and  computing  intercept  trajectories.  In  addition,  while  the 
velocity  of  the  target  is  assumed  constant,  the  velocity  of 
the  missile  is  assumed  to  vary.  Here,  the  tangential  accel¬ 
eration  is  assumed  a  positive  constant  while  thrusting  and 
a  negative  constant  while  coasting. 


As  a  first  step,  the  analysis  is  carried  out  in  two  di¬ 
mensions  with  the  hope  that  some  insight  in  the  three- 
dimensional  problem  will  be  achieved.  The  resulting  time- 
to-go  algorithm  is  tested  in  a  six-degree-of- freedom  simula¬ 
tion  by  projecting  the  current  missile  velocity  vector  onto 
the  plane  of  the  line  of  sight  and  the  target  velocity.  The 
optimal  intercept  time  is  computed  in  this  plane  and  used 
as  the  time-to-go  for  the  linear-quadratic  guidance  law. 

Optimal  Intercept  Problem 

Shown  in  Fig.  1  is  the  geometry  of  the  intercept  problem. 
The  XY  coordinate  system  represents  an  inertial  frame,  and 
the  X  axis  is  along  the  line  of  sight  at  t  =  0.  The  constant- 
velocity  target,  located  at  A’r  =  X,  at  t  =  0,  is  moving 
along  a  straight  line  which  makes  an  angle  0  with  respect 
to  the  X  axis.  The  missile  is  launched  at  an  angle  relative 
to  the  X  axis,  and  the  velocity  direction  0(t)  is  changed  by 
controlling  the  normal  acceleration  An(t).  U  x  =  Xt  —  Xu 
and  y  =  Yt—Yu,  the  equations  of  motion  of  the  engagement 
in  relative  coordinates  are  given  by 


X  =  Vt  cos  4  —  Vu  cos  0 
y  =  Vj-ain^  —  Vu  sin  B 
0  =  anIVst  ■ 


(1) 

(2) 

(3) 


The  tangential  acceleration  history  of  the  missile  is  as¬ 
sumed  to  be  constant  at^  >  0  while  thrusting,  that  is, 
for  t  <  te  where  is  fbe  known  engine  cutoff  time,  and  a 
constant  a,^  <  0  while  coasting.  As  a  consequence,  the 
velocity  of  the  missile  while  thrusting  is  given  by 


where 


kif  =  «i  +  ,  <  <  fe 


ui  =  Vm.  1  *1  =  <*(„ 


(4) 


(5) 


Similarly,  the  velocity  of  the  missile  during  the  coasting 
phase  is  given  by 


=  oj  -i-  hjt  , 


(6) 


where 


oi  =  Vm,  +  (ai«.  -  .  h  =  «(.<.  (7) 

The  prescribed  initial  conditions  are 

I  Vq  ”  Vog  ,  ~  (3) 


maximum  normal  acceleration  followed  by  zero  normal  ac¬ 
celeration,  that  is,  a  bang-bang  ronlrol.  The  bang-bang 
control  is  used  here  to  generate  a  minimum- time  trajectory, 
and  the  flight  time  is  used  as  a  prediction  of  the  time-to-go 
for  linear-quadratic  guidance  rules. 

In  the  solution  of  the  problem,  it  is  found  that  a  di¬ 
rect  intercept  can  be  achieved  for  Oza  <  <  Oua  (see 

Fig.  2)  where  Bza  is  the  initial  angle  for  the  zero-normal- 
acceleration  intercept  and  Qua  >s  the  initisJ  angle  for  the 
maximum-normal-acceleration  intercept.  For  >  Qua, 
the  missile  passes  in  front  of  the  target  during  the  maxi¬ 
mum  normal  acceleration  phase.  Then,  it  performs  a  360 
deg  turn  before  it  goes  for  the  intercept.  When  this  hap¬ 
pens  the  minimum  time  is  taken  as  the  time  to  the  point  of 
closest  approach.  A  similar  discussion  holds  for  Be  <  Bza- 

In  the  development  of  the  equations,  there  are  four 
important  times:  the  initial  time  to,  the  final  time  t/, 
the  engine  cutoff  time  U,  and  the  switch  time  t,  between 
“n  =  =  0-  It  is  assumed  that 


te<t.<  tj 


(10) 


so  that  o„  =  all  the  way  if  t,  =  ty  and  a,  =  0  all 
the  way  if  t,  =  te-  In  the  development  of  the  equations,  the 
engine  cutoff  time  is  assumed  to  satisfy  the  inequality 


f.  <fc<f/  • 


(11) 


Then,  if  te  <  (coasting  all  the  way),  the  correct  equations 
can  be  obtained  by  setting  te  =  and  if  >  </  (thrusting 
all  the  way),  tg  is  set  equal  to  (/. 

Constant  Normal  Acceleration 

The  equations  of  motion  (1)  through  (3)  can  be  inte¬ 
grated  for  the  case  where  is  constant.  These  solutions 
are  valid  for  the  cases  where  or  On  =  0. 

Since  the  missile  velocity  has  the  general  form  Vu  = 
o*  -I-  hst  where  k  =  1  for  thrusting  and  k  =  2  for  coasting, 
Eq.  (3)  can  be  integrated  as 

The  subscript  p  denotes  a  generic  starting  point;  it  could 
be  the  initial  point,  the  engine  cutoff  point,  or  the  switch 
point. 

Next,  with  Eq.  (12),  E)q.  (1)  can  be  integrated  to  yield 
x  =  x,  +  A(t,a„,k,t„0^)  -  A(te,a„,k,t„B,)  (13) 


where  the  subscript  s  denotes  a  specific  value.  Intercept  at 
the  final  point  requires  that 

xy  =  0  ,  yy  =  0  .  (9) 

If  Be  were  free,  the  control  for  the  minimum-time  tra¬ 
jectory  would  be  On  =  0,  that  is,  a  straight  line.  On  the 
other  hand,  if  Be  were  prescribed,  the  minimum-time  control 
would  be  infinite  normal  acceleration  to  rotate  the  velocity 
vector  instantaneously  to  the  above  straight-line  followed 
by  zero  normal  acceleration.  If  a  bound  were  applied  to 
the  normal  acceleration,  the  optimal  control  would  become 


where 

A(t,a„,k,t,,B,)  =  Vrcos^t-2ii^f^- 

•{cosK  +  !^^n(^)j 

Integration  by  parts  is  used  to  obtain  the  second  term  in 
Eq.  (14). 

Finally,  Eq.  (2)  for  y  is  integrated  in  the  same  way  as 
Eq.  ( 1 )  and  leads  to 

y  =  yp  +  f?(t,a„k  <„<;,)  -  f?(t„a„,k,t„tf,)  (15) 


where 


•{«“K  +  ff 'n(^)) 

+  (16) 

Zero-Normal-Acceleration  Intercept 


Og  >  OzAt  <»ii  ^  0-  Then,  the  equations  .pf  motion  are  ijite- 
grated  by  Euler  int^ation  for  a„  =  o,^  over  the  entire 
trajectory.  At  the  point  whw  xj  =  6,  the  value  of  y}’is 
Investigated  to  determine  whether  an  intercept  can  occur. 
If  the  missile  passes  behind  the  tuget  (y/  >  0)^  an  intercept 
can  be  achieved,  and  if  the  missile  passes  ahead  of  the  target 
[y/  <  0),  no  intercept  is  possible.  Should  the  missile  turn 
have  a'  small  raough  radius  that  2/  ^  0  after  a  reasonable 
time,  an  intercept  is  deemed  possible. 


In  this  section,  the  zero-normal-acceleration  intercept 
(straight-line  intercept  0  =  is  derived.  It  is  the  dividing 
line  between  the  a„  >  0  {0g  <  0za)  intercepts  and  a„  < 
0  (®o  >  ^za)  intercepts.  For  the  time  being,  the  cutoff 
time,  which  is  known,  is  assumed  to  satisfy  the  inequality 
constraint  t,  <  <  ty. 

If  a„  =  0,  Eqs,  (12),  (13)  and  (15)  can  be  applied  at 
and  if  and  the  results  combined  to  yield 


=  00  (17) 

xj  =  io-|-A(<„0,l,lo,<lo)->l(t„,0,l,t.,d.) 

-t-/l(<y,0,2,t„fl.)  -  A{tj,Q,2,to,0o)  (18) 

yj  =  y.  +  B(t.,0,l,l.,<?.)-B(l.,0,l,l.,<l.) 

-f-B(l/,0,2,l.,fl,)  -B(ty,0,2,t„tf.)  (19) 


For  intercept,  i/  =  y/  =  0  so  that  Eqs.  (17)  through  (19) 
involve  two  unknowns:  tf  and  0,  s  Oza-  These  equations 
can  be  solved  for  008  0^,4  and  sin  0^,4  as 


COS02A 


sin  02  A 


Vtom* 

VT«iwS  <>•*■¥« 


(20) 


which,  in  turn,  can  be  squared  and  added  to  obtain  the 
single  equation  for  ( / 


r  VyoK*  «>4t. _  , 

_ Vr««S  «>■>•»« _ *  _ 


(21) 


Once  t/  is  known,  0za  follows  from  Eqs.  (20).  If  to>tf, 
the  proper  equations  can  be  obtained  from  Eqs.  (20)  and 
(21)  by  setting  t*  =  f/.  On  the  other  hand,  for  <  <„  set 

tc  = 


Intercept  IVajectories 

For  an  intercept,  the  switch  time  is  assumed  to  satisfy 
the  inequality  tg  <  t,  <  tj  because  i.  =  ig  is  the  zero- 
normal-acceleration  intercept,  and  t,  =  tf  is  the  maximum- 
normal-acceleration  intercept.  There  are  four  possible  con¬ 
figurations  for  an  intercept  trajectory  depending  on  the 
value  of  the  known  engine  cutoff  time,  that  is,  U  <  tg, 
lo  ^  te  ^  t»i  <4  £  tfi  fc  >  1/.  The  equations  for 
U  S  tg  can  be  obtained  from  those  for  tg  <  tg  <  t,  by  set¬ 
ting  tg  =  Ig,  and  the  equations  for  tg  >  tf  can  be  obtained 
from  those  for  t,  <  tg  <  tf  by  setting  tg  =  tf.  Hence,  only 
two  sets  of  equations  need  to  be  derived. 


to<tc<t» 


For  this  case,  Eqs.  (12),  (13),  and  (15)  can  be  applied 
at  tg  and  t,  and  combined  to  yield 

+  ,  0.  =  0..hf;^fn(5f±hl^)  (22) 

X$  —  2»  "F  A(fef  Un,  1>  fp,  0o)  ~~  ^(foi  Un,  1,  lot  ^o) 

-FA(<.,o„2,fp,0e)  -/t(<p,o„,2,fp,0c)  (23) 

y»  ”  y*  "F  Itfpt^p)  ^(fpi  ^ni  1,  fpi  0») 

-FB(l.,o„,2,tc,0c)  -  B(fc,o„,2,t„0.)  (24) 

Along  the  straight-line  part  (a,  =  0),  Eqs.  (12),  (13),  and 
(15)  become 

0/  =  0.  (25) 

X,  =  X.  +  A{t,,O,2,t.,0,)-  A{t.,Q,2J.,0.)  (26) 

y/  =  y.  +  B(t/,O,l,f.,0.)-B(<.,O,2,f.,0.)  (27) 

Because  of  Eqs.  (14)  and  (16),  the  explicit  forms  of  Eqs. 
(26)  and  (27)  are  given  by 


Mazimum-Normal-Acceleration  Intercept 

The  maximum-normal-acceleration  intercept  is  the  di¬ 
viding  line  between  intercept  trajectories  and  closest-point- 
of-approach  trajectories.  It  is  possible  to  determine  this 
trajectory  by  applying  Eqs.  (12),  (13),  and  (15).  However, 
a  system  of  two  equations  in  two  unknowns  results,  and 
while  it  can  be  solved  efficiently  with  careful  coding,  a  sim¬ 
pler  approach  is  used  to  determine  whether  an  intercept  is 
possible. 

Once  OzA  is  known,  the  actual  0,  can  be  compared  _witji 
OzA  to  determine  the  sign'd!  a„:  0,  <  Oza,  >  0  or 


1/  =  I.  -F  {VrCDsd  -  (oj  -F  56,(1/  -  l,)]cos0,}(l/  -  1.X28) 
yj  =  y.  +  {Vrsind-  K  +  ^6,(1/  -l,)]sin0,}(l/  -  1,)(29) 

The  solution  process  is  to  vary  1,  until  x /  and  y^  are 
both  zero.  This  is  accomplishet^y  setting  z/  =_0  _and 
solving  Eq.  (28)  analytically  for  tf.  Eq.  (27)  is  a  quadratic 
equation  and  has  two  roots.  The  correct  root  is  the  smallest 
value  of  l/'which  is  Larger.than  1..  Then,  ^hejolution  for 
tf  is  substituted  into  Eq.  (29)  where,  because  t,  is  guessed, 
Hf  ^  0.  Hence,  1,  is  varied  until  y/  =  0.  In  this  way,  it  is 
only  necessary  to  solve  one  equation  in  one  unknown. 


ta<tc<tf 

If  the  engine  cutoff  time  is  between  the  switch  time  and 
the  final  time,  Elqs.  (12),  (13),  and  (15)  lead  to 


0,  =  0.,  0.=  00  + 


6i  \aj  +  5il«/ 


Zc  —  Z®  +  A(f«,  flni  li  foi  ^o)  ^(f«i  1»  fo» 

+A{U,Q,l,t.,0.)  -  A{i.,0,l,t.,0.)  (31) 

ye~  Vo  +  ^(fi»  ®i»i  li  fo»^«)  ^(fdi  ®n»  1»  f**f^o) 

+B{u,  0,  Sit.,  0. 1,  t..  0.)  (32) 

For  the  coast,  the  corresponding  equations  are 

01  =  0c  (33) 

x,  =  Xo  +  A{tj,  0, 2,  U,  0.)  -  AiU,  0, 2, 0c)  (34) 

yf  =  yc  +  Bitf,Q,2,to,0o)  -  Bito,O,2,to,0o)  (35) 

The  explicit  forms  of  Eqs.  (34)  and  (35)  are  similar 
to  those  of  Eqs.  (28)  and  (29).  Hence,  the  same  solution 
process  can  be  used. 

Closest-Point-of-Approach  Trajectories 

For  a„  =  an«„  over  the  entire  trajectory  and  for  t,  < 
tc  S:  tf,  Eqs.  (12),  (13),  and  (15)  can  be  applied  at  t.  and 
t /  and  combined  to  obtain 

0/  =  ec+i^in{^)  ,  =  +  (36) 

Zy  —  Xo  +  A(fei  0„,  1,  t.,  0.)  —  i4(ttf,  On*  1*  tc,0o) 

+Aitf,  2,  to,  0.)  -  Aitc,  a,,  2,  0^)  (37) 

y/  “  y«  +  Bitf,  a„,  1,  fg,  0.)  —  BHo,  On,  1»  to,  0o) 

+Bit,,a„,2,to,0o)  -  Bito,ao,2,to,0o)  (38) 

The  equations  for  the  case  where  t.  <  t.  are  obtained  from 
Eqs.  (36)  through  (38)  by  setting  te  =  t,.  Similarly,  for 
t.  >  tj,  set  tc  =  tj. 

The  point  of  closed  approach  is  obtained  by  minimizing 
the  performance  index 

d*  =  z*  +  y*  .  (39) 

This  is  accomplished  by  solving  the  algebraic  equation 


by  bisection.  The  derivative  in  (40)  can  be  taken  analyti¬ 
cally.  To  verify  that  the  solution  is  a  minimum,  the  second 
derivative  d^<P/dt^j  is  checked. 

Numerical  Results 

The  algorithm  followed  in  computing  a  trajectory  is  the 
following: 

1.  Given  Vr,  <(>,  Vmo,  a,„..,  a,^,  t.,  z„  y.,  and  t.,  com¬ 
pute  0ZA-  This  determines  whether  a,  >  0  or  a,  <  0. 


2.  Given  <in„„  and  0.,  determine  whether  or  not  an  in¬ 
tercept  can  occur. 

3.  If  an  intercept  is  possible,  calculate  the  switch  time 
for  z/  =  y/  =  0. 

4.  If  an  intercept  is  not  possible,  calculate  the  final  time 
for  closest  point  of  approach. 

Some  results  have  been  obtained  for  the  following  condi¬ 
tions: 

Vt  =  500  ft/sec,  0  <  <  ISO  deg,  Vm.  =  1,000  ft/sec 

=  25  g's,  =  -10  g's,  to  =  0  see, 

Xo  =  5, 000  ft,  yo  =  0  ft,  tc  =  2.6  sec, 

“n...  =  100  g's,  -70  <  <?.  <  70  deg  (41) 

which  are  typical  of  a  conceptual  bank-to-turn  missile  en¬ 
gagement.  The  limit  on  0.  is  representative  of  the  field  of 
view  of  a  passive  seeker.  TVajectories  are  presented  in  Figs. 
3  through  7  for  different  values  of  ^  primarily  to  illustrate 
that  the  computation  process  is  valid. 

The  time-to-go  algorithm  has  been  tested  in  a  six- 
degree-of-freedom  simulation  of  a  bank-to-turn  missile.  The 
missile  velocity  vector  is  projected  onto  the  plane  of  the  line 
of  sight  vector  and  the  target  velocity  vector;  the  time-to-go 
is  calculated  in  this  plane;  and  the  result  used  as  the  time- 
to-go  for  three-dimensional  flight.  Miss  distance  results  are 
shown  in  Table  1  for  a  number  of  engagement  senerios.  The 
results  labeled  Case  I  are  for  the  case  where  time-to-go  is 
calculated  as  range  divided  by  closing  speed.  The  results 
labeled  Case  II  are  for  the  case  where  bang-bang  trajecto¬ 
ries  are  used  to  calculate  time-to-go.  In  six  engagements, 
the  bang-bang  formula  allows  the  missile  to  hit  the  target 
(within  10  ft)  when  the  range-over-closing-speed  formula 
allowed  a  miss.  Only  once  does  the  reverse  happen.  In  all 
other  engagements  both  formulas  generate  a  hit  or  a  miss 
simultaneously.  Overall,  the  bang-bang  formula  (C^Lse  II) 
provides  better  results. 

Discussion  and  Conclusions 

A  method  for  predicting  time-to-go  for  homing  missiles 
using  a  linear-quadratic  guidance  rule  has  been  developed. 
It  is  based  on  minimum-time  trajectories  for  senerios  when 
intercepts  are  possible  and  on  minimum-miss-distance  tra¬ 
jectories  when  intercept  trajectories  are  not  possible  in  a 
reasonable  time.  The  missile  velocity  is  variable  in  that 
the  tangential  acceleration  has  a  known  positive  value  when 
thrusting  and  a  known  negative  value  when  coasting  and  the 
engine  cutoff  time  is  known.  On  the  other  hand,  the  target 
velocity  is  assumed  constant. 

In  general,  the  optimal  trajectories  are  composed  of  con¬ 
stant  normal  acceleration  segments.  For  an  intercept  trajec¬ 
tory  the  control  is  maximum  normal  acceleration  followed 
by  zero  normal  acceleration.  For  a  minimum-miss  trajec¬ 
tory,  the  control  is  maximum  normal  acceleration  all  the 
way.  Numerical  results  for  the  time-to-go  algorithm  have 
been  presented  to  verify  the  computational  procedure.  Also, 


the  time-to-go  algorithm  has  been  tested  in  a  six-degree-of- 
freedom  simulation.  In  general,  the  proposed  time-to-go  al¬ 
gorithm  produces  better  miss  distances  than  those  obtained 
from  the  range-over-closing-speed  method. 
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Fig.  1  Intercept  geometry 


Table  1:  Simulation  Results 
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Fig.  2  Typical  Trajectories 
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Fig.  3  Trajectories  for  =  0  deg 
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'ig.  4  Trajectories  for  ^  =  45  deg 
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Fig.  7  Trajectories  for  ^  =  180  deg 
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'ig.  5  Trajectories  for  ^  =  90  deg 
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